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Abstract 

The problem of measuring the resistivity of isotropic 
samples of finite dimensions in the form of rectan¬ 
gular parallelepipeds using the four-probe technique 
was considered. Two variants of contact arrange¬ 
ments were studied: (1) four collinear probes are po¬ 
sitioned on one side of a sample symmetrically with 
respect to the other sides, and (2) two probes on one 
side of a sample and two on the opposite side are 
placed precisely in opposite positions and symmet¬ 
rically with respect to the other sides of the sample 
(the Schnabel method). Solutions of the problem of 
the electric field potential distribution in a sample for 
different positions of the current contacts were found. 
The solutions were obtained in the form of double se¬ 
ries and methods of their summation are presented. 
The obtained results are extended to the case of mea¬ 
suring the resistivity of anisotropic samples when the 
resistivity tensor has two independent components. 
The results of using the developed technique for mea¬ 
suring the resistivity of such a highly anisotropic ma¬ 
terial as highly oriented pyrolitic graphite using the 
Schnabel method are presented. 


1 Introduction 

Four-probe methods are widely used to measure the 
resistivity of materials. The essence of these methods 
is as follows. Four point contacts are placed on the 
surface of an investigated sample, a preset current 
/ from an external source is passed through two of 
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them (current contacts), and the potential difference 
AV is measured between the two other (potential) 
contacts. It is obvious that the measured potential 
difference is proportional to the current value /: 

AV = RI , (1) 

where the coefficient R, which will be called the con¬ 
ditional resistance, is a complex function of the sam¬ 
ple geometry, the positions of the current and poten¬ 
tial contacts on it, and, certainly, the resistivity p of 
the sample material. 

If the studied material is isotropic and is avail¬ 
able in the form of long wires or thin long strips, the 
problem is trivial. The current contacts are placed 
at the ends of the sample, and the current is then 
distributed uniformly over the sample cross section 
S' at a sufficient distance from the sample ends. The 
potential contacts are mounted at a rather long dis¬ 
tance from the current contacts at a specified dis¬ 
tance l from each other, and the potential difference 
AV between them is measured at a given current I. 
The resistivity p of the material is determined from 
the formula R = A V/I = pl/S. This technique is 
well known and widely used in practice. 

In more complex cases where the investigated ma¬ 
terial is represented by, e.g., small crystals (experi¬ 
mental solid state physics), massive semiconductor 
silicon and germanium wafers (industrial produc¬ 
tion in microelectronics), massive alloy blocks (met¬ 
allurgy), etc., the problem becomes much more diffi¬ 
cult. To determine p from the measurement results, 
it is necessary to theoretically calculate the form of 
the function R in formula CD for a particular ge¬ 
ometry of the experiment. This is a rather complex 
problem, especially if the resistivity p is anisotropic. 
However, this problem can be successfully solved for 
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samples with a definite shape at a special arrange- 2 
rnent of the measuring electrodes on the sample sur¬ 
face. 


Isotropic Sample of Finite Di¬ 
mensions 


The key moment is the solution of the Laplace 
equation for the electric field potential u(x, y, z ) 
in the investigated material with the correspond¬ 
ing boundary conditions. In the case of an infinite 
isotropic conducting half-space or an infinite plate, 
the corresponding solutions are known well and can 
be found, e.g., in [Mi¬ 
lt usually happens in practice that the studied ob¬ 
jects are samples of small dimensions in which the 
distance between the measuring contacts is compa¬ 
rable to the sample size. In this case, the results 
described in m are inapplicable, and a solution of 
the Laplace equation for a sample with finite dimen¬ 
sions should be sought. This problem can be solved 
analytically for samples that are shaped as a rectan¬ 
gular parallelepiped |5H2|. This solution will be pre¬ 
sented below in the form that is, in our opinion, more 
convenient for application than that in [5H3- The 
matter is that the solutions that are obtained using 
the method of separation of variables have the form 
of double series, which, as a rule, poorly converge. 
It is shown below how such series must be correctly 
summed and correct results can be obtained. The 
technique we proposed in m allows consideration 
of samples with finite values of all three dimensions. 

In the present paper, we make an attempt to ana¬ 
lyze in detail the theoretical questions of such mea¬ 
surements for isotropic substances with the subse¬ 
quent application of the obtained results to a cer¬ 
tain class of anisotropic materials. In connection 
with this, it should be noted that, e.g., in solid state 
physics, exactly highly anisotropic systems with lay¬ 
ered structures, such as high-temperature supercon¬ 
ductors, bismuth-chalcogenide-based topological in¬ 
sulators, highly oriented pyrolitic graphite, etc., at¬ 
tract the greatest interest. The proposed technique 
with the use of four point contacts is especially con¬ 
venient for such materials and allows one to obtain 
the absolute values of the longitudinal (along layers) 
and transverse (across layers) resistivities and, cor¬ 
respondingly, the value of the resistivity anisotropy. 

Two variants of the contact arrangement are con¬ 
sidered: (1) four collinear probes are positioned on 
one side of a sample, and (2) two contacts on one 
side of a sample and two on the opposite side are 
placed precisely in opposite positions (the so-called 
Schnabel geometry mm- 


Let us consider an isotropic sample in the form of a 
rectangular parallelepiped with the sides a and b and 
thickness d\ the directions of the coordinate axes are 
shown in Fig. [I] 


a 



Figure 1: Sample of finite dimensions in the form of a 
rectangular parallelepiped with the length a, width 6, and 
thickness d and the corresponding directions of the x, y , 
and z coordinate axes (an isotropic sample) and x\, X 2 , 
and X 3 axes (an anisotropic sample). The coordinate ori¬ 
gin is at the center of the parallelepiped. 


2.1 Two Current Contacts Symmetrically 
Positioned on One Side of a Sample 


Let us first consider the case where point current 
contacts are placed on one side of the sample: z = 
—d/2. To simplify this problem, we begin to consider 
a practical case where both current electrodes are 
symmetrically positioned on this side (Fig. 2a); i.e., 
their coordinates are (±xq,0, —d/ 2 ). 

If a stationary current / flows through the contacts 
A, B, an electric field establishes in the sample. Its 
potential u(x, y, z) satisfies the Laplace equation 


. d 2 u d 2 u d 2 u 

u{x, y, z) = "x o + + yr~2 = 0 

dx z dy z dz z 

with the boundary conditions 


( 2 ) 


du 

dx 

du 

dz 


du 
= °, w- 
x=±a /2 d V 


y=±b /2 


du 
= °’ ITz 


z=d/2 


= 0, 


z=—d/ 2 


= -pl[5(x + x 0 ) - 5(x - x 0 )]<%). 


(3) 

The first three conditions in (J3J) designate the ab¬ 
sence of a current through the interfaces between the 
sample and environment, and the last condition cor¬ 
responds to the point nature of the contacts through 
which the measurement current / flows. 
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(a) 



(b) 


(c) 


(d) 


Figure 2: Geometry of the arrangements of the contacts 
on the sides of a sample in the form of a rectangular paral¬ 
lelepiped in the four-probe method: (a) two current con¬ 
tacts A, B with the arriving and outgoing current I, re¬ 
spectively, are positioned symmetrically on one side of the 
sample, and if the potential contacts D , C are positioned 
symmetrically with respect to the current contacts on the 
opposite side, we have the so-called second Schnabel ge¬ 
ometry; (b) four current contacts with the entering and 
outgoing current I /2 (as shown with arrows) are posi¬ 
tioned on opposite sides of the sample; (c) four current 
contacts with the entering and outgoing current 1/2 (as 
shown with arrows) are positioned on one side of the sam¬ 
ple; and (d) the first Schnabel geometry (the current and 
potential contacts are, respectively, at the points A, D 
and B , C). 


The standard method of separation of variables [3J 
consists in the fact that the solutions of Eq. (0) are 
sought in the form of the product of three functions 
u(x,y,z ) = X (x)Y (y) Z (z), each of which depends 
only on one variable. The functions X(x), Y(y), and 
Z(z) satisfy second-order ordinary differential equa¬ 
tions 


d?W(w) 
dw 2 


— cW (w) = 0 . 


which, depending on the value and sign of the con¬ 
stant c, have one of three fundamental solutions: 
W(w) = C\ w + C‘ 2 : W(w) = A\ sin A w + B\ cos A w, 
or W(w) = Die Xw + D 2 e~ Xw . 

The general solution of Laplace equation d^J) is 
composed of partial solutions of this kind. Boundary 


conditions Q and the symmetry considerations al¬ 
low the number of acceptable partial solutions to be 
significantly reduced. It was shown in m that for 
the given geometry, the solution of Laplace equation 
§2\) with boundary conditions ([ 2 ]) is described by the 
expression 


u(x,y,z) 

8 pld cosh 7(2 — d/ 2 ) . (2k + l)vrxo 

- \ -sm- 

ab ]~/ 0 ydsinhyd a 


(2k + l)nx 2mry 

x sm-cos —-— , 

a b 

where the quantity 7 is defined by the formula 


(4) 


7 kn — Tt 



2 


+ 



(5) 


and the coefficient 6 n is 


f 1 , n^O, 

\ 1 / 2 , n = 0 . 


( 6 ) 


This result can be immediately applied to mea¬ 
surements of the resistivity p both in the second 
Schnabel geometry (Fig. [2Jr) and in the first collinear 
geometry, when the outer (A and D ) and inner (B 
and C) contacts serve as the current and potential 
contacts, respectively, in a similar way, as is shown in 
Fig. for another case. However, this will be done 
later, and let us now consider two auxiliary problems 
with four current contacts. 


2.2 Four Current Contacts on Opposite 
Sides of the Sample 


Let us first consider the auxiliary problem (Fig. [2}o) 
where there are four symmetrically positioned cur¬ 
rent contacts; thus, the contacts A and B have the 
coordinates (±xo, 0, — d/2), and the contacts C and 
D are positioned at the points (±xo, 0, d/2). The 
current 1/2 enters each of the contacts (A. B), and 
the same current 1/2 goes out of the opposite con¬ 
tacts (C, D). 

We must solve Laplace equation © with the 
boundary conditions 


du 
dx 
du 
dz 


du 

= 0 ’ 7T 

x=±a /2 °y 


= 0, 


y=±b/2 


Pl< 


z=±d/2 


= -—[K x + x o) + H x 


XO )]<%). 


From the symmetry considerations, the solution 
for the potential u(x, y , z) must be symmetrical with 
respect to the arguments x and y. 
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In this case, the solution of the Laplace equation 
for the potential u(x, y, z) has the form [10] 


u(x, y, z) 

pi Apld OkOnSmh^z 

ab 7dcosh(7<i/2) 


ab 


2 knxo 2kirx 2niry 

x cos- cos-cos —-— 

a a b 

where 7 is defined by the formula 


7 kn = VTA 


Zk\ 2 f 2n\ 2 


a ) 


b 7 ’ 


( 7 ) 


( 8 ) 



and 9 k and 9 n are defined by formula ((6]) . The prime 
in the summation sign in formula © means that the 
term with k = n = 0 is excluded from the summa¬ 
tion. 


2.3 Four Current Contacts on One Side 
of the Sample 


Let us consider one more auxiliary problem, when 
there are four current contacts that are positioned 
symmetrically along a straight line on one side of 
the sample, as is shown in Fig. [ 2 ]x The contacts 
B,C and A, D have the coordinates (±21, 0 , — d/ 2 ) 
and (±22, 0 , — d/ 2 ), respectively. The current //2 
enters each of the contacts A and D, and the same 
current leaves the contacts B and C. 

In this case, we must solve Laplace equation (© 
with the boundary conditions 




du 

dx 

du 

dz 


du 
= °, tt 

x=±a /2 d V 


V=±b /2 


_ du 
= °’ ~d~z 


= 0 , 


z=d/2 


z=-d/2 


[^( x + £2) — S(x + xi)— 
— 5 (x - 21) + S(x - x 2 )]d(y). 


Assuming for definiteness that the contacts are po¬ 
sitioned equidistantly, i.e., 21 = s /2 and 22 = 3 s/ 2 , 
we obtain the solution for the potential distribution 
in this problem in the form m 


u(x,y,z) 

8 pld 9k9 n cosh 7(2 — d/2) . kirs . . 

=-w- / -;-—rrx- sin- ( 9 ) 

ab f ~—/ 0 7dsmh(7d/2) a 

2k-KS 2/c7T2 2mry 

x sin- cos - cos —-— , 

a a b 

where 7 is defined by formula ©, and and 6 k and 
9 n , by formula d6j). 

Now, the obtained results can be applied to prac¬ 
tical measuring schemes. 



Figure 3 : (a) Summation of the double series along the 
diagonals on the (k, n) plane; (b) the result of calculating 
the function H as a function of the considered diagonals 
for a sample with dimensions of a = 9 mm, 6 = 3 mm, d 
= 0.3 mm, and s = 3 mm. 


2.4 Measurements using the Schnabel 
Method. Problem of Series Summa¬ 
tion 

Let us first consider the simpler second Schnabel 
geometry (Fig. ©) in which the current contacts 
have the coordinates (±s/2,0, —d/2) and the poten¬ 
tial contacts are placed on the opposite side at the 
points (±s/2,0, d/2) symmetrically with respect to 
the current contacts. 

Substituting these values of the arguments into for¬ 
mula ( 4 ), we easily find the potential difference Voc 
between the potential contacts: 


Vdc 


* n T 1 00 n -2 (2fc+lY 

16 pld On sin — ^~ 

ab s i n h 7 ^ 


( 10 ) 


Thus, the conditional resistance R2 can be ex- 


4 

































pressed by the formula 


7 


7/2 7/2 


R2 


Vdc 

I 




( 11 ) 


where the function H is defined by the expression 


H = 


16 d 2 y, Q nS in 2 ®+g™ 
ab pd sinh z>d 

k,n =0 


( 12 ) 


in which the quantity 7 kn is found using formula © • 

Hence, the resistivity p for an isotropic sample can 
be found from formula ED using the results of one 
measurement of the resistance R2 in the second Schn¬ 
abel geometry. As is seen from formula ([ 1211 . the 
function H depends only on the sample dimensions 
and the distance s between the probes and can be 
easily calculated using this formula via summation of 
the double series. This series perfectly converges ow¬ 
ing to the presence of exponents in the denominator 
of its common term (hyperbolic sine). Summing on 
the plane of the coefficients {k, n) can be performed 
by different methods: in squares (k < N,n < N ) or 
in circles (k 2 +n 2 < N 2 ), where the chosen size of the 
summed area (number N) is rather large. However, 
the most suitable method is to sum along diagonals 
(k + n = N ), as shown in Fig. ( 3 jr. The summation 
procedure terminates at a sufficiently large N, when 
the required accuracy is attained. As a rule, it is 
sufficient to sum terms on the first several tens of 
diagonals on the (k, n) plane. Figure [ 3 }o shows, as 
an example, the change in the result of calculating 
the function H as a function of the number of con¬ 
sidered diagonals N for a sample with dimensions of 
a = 9 mm, 6 = 3 mm, d = 0.3 mm, and s = 3 mm. 

Let us now consider the problem of the elec¬ 
tric held potential distribution for the arrange¬ 
ment of the current and potential probes in accor¬ 
dance with the first Schnabel geometry (Fig. [ 2 ji) in 
which the current contacts are placed at the points 
(—xo, 0, ±d/2), and the potential contacts are at the 
points (xo, 0 ,±d/ 2 ). In this case, the sought distri¬ 
bution can be obtained via superposition of three 
fields in accordance with Fig. [ 4 l Namely, the desired 
potential can be represented as the sum u(x, y, z) = 
u±(x, y, z) + U2(x, y, z) + us(x, y, z), where the poten¬ 
tial u\ corresponds to the problem with four current 
contacts (considered in Subsection 2 . 2 ), which are 
positioned at the same points (Fig. [ 2 }o) . The poten¬ 
tial U2 corresponds to the second Schnabel geometry 
(Fig. [ 2 ji) with the measuring current 1 / 2 , and the po¬ 
tential 1x3 corresponds to the same second Schnabel 
geometry with the current I/2 and current electrodes 
positioned on the side z = d/2. 


Ill 7/2 



7 IH 7/2 III 7/2 


Figure 4 : Potential in the first Schnabel geometry as a 
superposition of three known fields. 


The potential u\ is expressed by formula ©, and 
the potential U2 is expressed by formula @ in which 
1/2 must replace the current value I. The potential 
113 is expressed by the same formula as U2 in which 
the signs of x and z are reversed, i.e., uz{x,y,z) = 
U2{—x,y,—z). As a result, the following expression 
for the potential in the sample is obtained: 

u{x,y,z) 


ab Z 




ab k % o yidcosh^ 




a 


x cos 


2k-KX 2mry Apld 

- c °s —-— + —— 

a b ab 


x 


x cos 


x cos 


E °° sinh72Z ( 2 A: + l)7rxo (2fc + l)7r.T 

-7 sin- sm- 

t,n=o 72 d cosh ° a 

2niry Apld 

h nh 


k,n =0 _ 


smli7iz 2kirxo 

cos- 

a 


7idcosh =r 

2knx 2mry 0 n sinh722: . (2k + l) 7 rxo 

-cos —--—r sm- 

b 7 2 d cosh 2|- a 


a 


x sin 


(2k + 1 ) 7 tx 2nny 
-cos- 


a 


( 13 ) 


keeping in mind that for k = n = 0, the first term in 
the square brackets must be replaced by z/( 4 d) via 
a passage to the limit 71 —» 0 . Here, 



Now, we can calculate the potential difference Vbc 
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measured in the first Schnabel geometry between the 
potential contacts B and C (Fig. EJl). Assuming 
xq = s/2 (s is the distance between the contacts), 
x = s/2, and y = 0, we obtain from (fT3l) : 


Vbc =^yi 


ab 


k.n=0 L 


0k0 n , ,71 d 2 k7TS 

-- tanh- cos - 

7 ia 2 a 


On , ,72^-2 ( 2fc + 1 ) 7rs 

-- tanh- sin - 

7 2 d 2 2 a 


Thus, the conditional resistance R\ measured in 
the first Schnabel geometry is determined as 


d V BC P n 
Rl = — = d a 


(14) 


where the function G has the form 


G = 


8 d 2 
ab 


E 

k,n =0 


@k@n , i 2 

-- tanh- cos 

. 7i d 2 


kns 


--tanh 

7 2 d 2 


72 d . 2 ( 2 A; + l)7rs 


(15) 


sin 


2 a 


Unfortunately, it is impossible to use formula (1T51) 
because the double series in it converges poorly, and 
for the potential contacts in this geometry, which are 
positioned at the points (s/2,0, ±d/2 ) on the sample 
surface, the summation result does not correspond 
to the physical meaning of the problem. This is due 
to the fact that, in this case, the observation points 
correspond to places where the fictitious current con¬ 
tacts are located in accordance with Fig. [4j 

In order to get out of this difficult situation, we 
shall consider hypothetic potential contacts that, in 
contrast to real contacts, are located not at the sam¬ 
ple surface but inside it at the points (s/2, 0, ±(d/2 — 
e)), where e/d is a rather small value. The result we 
need for the contacts that are on the sample surface 
can be obtained at e tending to zero. 

Adopting this approach, we obtain the following 
expression for the measured potential difference Vbc 
from formula (1131) : 


Vbc = 


8 pld 


OO 


ab 

On 

lid 


k,n=0 


6 k 0 n sinhyid — e) 2 kns 


7i d 

sinli 7 2 (f - e) 


cosh ^ 


■ cos 


cosh ^ 


■ sm 


2 (2 k + 1)7TS 


2 a 


From here, as was above, we obtain the expression 
for the function G , which is written in a more suit¬ 
able form for practical calculations, by separating the 


term corresponding to k = n = 0 and expressing the 
hyperbolic sine and cosine through the exponents: 


8 nd 2 f 1 e a sinh^(| — e) . 2 7rs 

ab 1 8 4 d 2nd cosh ^ 2a 


OO f 


+ E 

k,n=0 


e k e n e~ 11£ - e-'lVd-e) kns 

-r-:--j- cos - 

7 i d 1 + e i ld a 


9 n e~^ £ - e~^ {d - £ ) (2k + lWs I 

-;- j - sin - > . 

72 d 1 + e 7 2d 2 a j 

(16) 


1.769 


1.768 


0.005 


0.01 

Z/d 


0.015 


0.02 


Figure 5: Dependence of the function G on the depth of 
embedding £ of the potential contacts. Numbers near the 
curve denote the approximate number of diagonals N on 
the (fc, n) plane. 


The double series in (1161) converges because of the 
presence of the corresponding exponents, and the 
convergence rate increases with an increase in the 
deepening e. However, to obtain the sufficient ac¬ 
curacy, one has to consider several hundred or even 
thousand diagonals on the (k. n) plane. 

As an example, Fig. [5] shows the plot that presents 
the results of calculating the function G from formula 
(1161) for a sample with dimensions of a = 9 mm, b = 
d = s = 3 mm at different e values. The approximate 
number of diagonals N on the plane of indices (fc, n), 
over which double series (1161) must be summed for the 
required accuracy to be obtained, is also indicated in 
this curve. This curve is a parabola from which it is 
easy to find the value of G at e = 0 by extrapolation. 
If two results G'i(ei) and G^^) are taken and the 
curve G{s) is approximated by the parabola G(e) = 
G( 0) + ge 2 , then we obtain the following value for 
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G(0): 

0,0) = Cl g-g e l (17) 

e 2 £ 1 

If £2 = 2ei, formula (ED transforms into G = 
(4G] — G 2 )/ 3 , and the result very weakly depends 
on £\. Just in this way, using two G values, the func¬ 
tion G( 0 ) was further calculated. 


2.5 Four Collinear Contacts on One Side 
of the Sample 

Let us now consider a widely used practical case 
where all the four probes are positioned on one side 
of a sample equidistantly and symmetrically along 
one line (e.g., as is shown in Fig. [2f). We define 
the first collinear geometry as the configuration of 
probes when the current contacts are at the points 
(A, D) and have the coordinates (±3s/2,0, —d/2), 
and the potential contacts ( B,C ) have the coordi¬ 
nates (±s/ 2 , 0 , —d/ 2 ), where s is the distance be¬ 
tween the neighboring probes. 

Using formula ©, we find the potential difference 
Vbc and the conditional resistance 


Ri 


Vbc 

I 


16 pd 
ab 

Vk 


E 


9 n cosh yd 
yd sinhyd 


k,n =0 

(2 k + 1)7TS 
2 a 


sin</y sin3</y, 


Thus, by analogy to formulas ED and ED, Ri 
can be written in the form 


where 


Ri = -M 

a 


16d 2 0 n cosh yd 

ab 7 d sinh yd 


sinyy sin3yy, 


(18) 


and the parameter y k n is determined by formula © . 
By measuring R\ and calculating the value of the 
function L\ using formula (11811 . we can find the ma¬ 
terial resistivity p. 

When the function L\ is calculated, we also en¬ 
counter a very poor convergence of the double se¬ 
ries in formula ED- When summing in squares 
(k < N,n < N), a partial sum of this series oscil¬ 
lates around a certain average value without damp¬ 
ing with an increase in N. To obtain a correct result, 
the summation should be performed along diagonals 
{k + n = N ) on the ( k , n) plane. In this case, these 
oscillations gradually decay with an increase in N 
and approach the average value. 


i 1 112 m in 



i in m in 



(b) 


II I 



Figure 6: (a) The second collinear geometry as a su¬ 
perposition of three known fields; (b) the third collinear 
geometry as a superposition of three known fields; and 
(c) the first collinear geometry as a superposition of the 
second and third geometries. 


To improve the convergence, we can use the same 
technique as that used in Subsection 2.4, i.e., con¬ 
sider that the potential probes are embedded into 
the sample to the depth e and have the coordinates 
(±s/ 2 , 0 , —d/2 + £); after that, we determine the de¬ 
sired value of the function L\, when e tends to zero. 
As a result, we obtain 


T 16 d 2 

Ll -~ 


a cosh^(d —e) 

M sinh — sin »’°‘ iin3 ^ 

a 


+ E 

k,n=0 


0 n e~^ £ + e-^d-A 


yi d 


l _ e -27i d 


sinyy. sin3</y. 


(19) 


where yi and </y are defined by the formulas 



(2k + 1)7TS 


Let us now consider the second collinear geom¬ 
etry, when a current is fed to the contacts A 
and C with the coordinates (—3s/2, 0, — d/2) and 
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(+s/ 2 , 0 , —d/ 2 ), and the potential difference is mea¬ 
sured between the points B and D with the coordi¬ 
nates (—s/2,0, —d/2) and (+3s/2,0, — d/2). In this 
case, the desired distribution of the electric field po¬ 
tential in the sample can be obtained as a superposi¬ 
tion of three fields in accordance with Fig. 6 a, i.e., by 
the same method as that applied in Subsection 2.4. 
The field potential u\ is described by formula Q, the 
potential 112 is given by formula @ at xq = 3s/2, and 
the potential U 3 is also described by formula J5J) but 
at xq = s/ 2 . 


the form 

16d 2 


L 3 — 


ab 
00 1 r 


a cosh f(d-e) . 2 2 

-^-sm po cos 2 ipo 


27t d sinh — 


a 


+ E 


k,n =0 _ 


9 k e n e-^ + e~^ 2d - £ ) . 2 , . 2n , 

-^d - 1-e-W - Sm ^ fcSm 2 ^ k 


On e-T' 1£ + e-T' 1 ( 2d - £ ) , 2 2 ' 

■ sin“ pk cos 2 pk 


71 d 1 - e~ 2 ^ d 


(23) 


After appropriate calculations are performed, the 
measured potential difference between the contacts 
B and D and the conditional resistance R 2 = 
Vbd/Iac ca n be calculated. The latter is expressed 
via the formula R 2 = (p/d)L 2 , where the function 
L 2 is expressed by the formula that is analogous to 

m 


T 16 d 2 

2 ab 


a cosh 5 (d — s) 9 9 

2 ^ si „ h ^ cos msin 2w 

a 


00 1 


-E 


OkOn e~ 72£ + e~ 72 ^ 2d ~ £ ^ 
72 d 1 — e~ 2 ^ 2d 


sin 2 -i/fc sin 2 2 i/k 


0 n e _71£ + e - 7 i( 2d - £ ) 
7 ~/d 1 - e ~ 2 ^ d 


n . o 

cos^ pk sin 2 pk 


( 21 ) 

where 71 and pk are defined by formulas ( 1201 ) . and 
the parameters 72 and ipf. have the form 


72 = 7T 


/ 



2 


+ 



V’fc 


kns 

a 


( 22 ) 


The third collinear geometry, where the current 
contacts are the points ( A,B) and the potential dif¬ 
ference is measured between the contacts (C, D), can 
be considered in a similar manner. The sought field 
distribution is found as a superposition of three fields 
in accordance with Fig. [SJo. The potential u\ is de¬ 
scribed by formula ©, the potential U 2 is given by 
formula (JU) at xq = 3s/2 and the current 1 / 2 , and 
the potential U 3 is also described by formula (JU) but 
at xq = s /2 and the current — 1 / 2 . 

Performing the same calculations, we obtain that 
the resistance R 3 measured in this geometry is ex¬ 
pressed by the formula R 3 = {p/d)L 3 , where L 3 has 


Here, 71 , 72 , Pk, and ipk are defined by formulas 
(1201) and (122b . 

The functions L\, L 2 , and L 3 are not independent 
and are interrelated via the expression L\ = L 2 + L 3 , 
which becomes obvious from Fig. [Bfc in which the 
first collinear geometry is presented as the superpo¬ 
sition of the third and second collinear geometries. 
This can be also demonstrated directly from formu¬ 
lae (USD, (12Tb . and (1231) . 

Thus, all the above-described four-probe tech¬ 
niques (collinear arrangement of the probes and the 
Schnabel method) allow one to obtain the resistiv¬ 
ity of an isotropic sample, which has the shape of a 
rectangular parallelepiped of finite dimensions, from 
one measurement 0 

3 Measurement of the Resistivity 
of Anisotropic Samples 

Anisotropic conducting substances are of great in¬ 
terest for solid-state physics. It was shown in m 
that the problem of the potential distribution in an 

1 In all previous sections, point contacts were considered. 
It is clear that from the experimental standpoint, this im¬ 
poses appreciable limitations on the current value that passes 
through the current contacts and, hence, reduces the sensitiv¬ 
ity of the method, because the measured potential difference 
between the potential contacts may be so low that it cannot be 
measured with the required accuracy. I 11 the study by Schn¬ 
abel [2], an attempt was made to extend the theory of the 
four-probe technique for measuring the material resistivity to 
the case of nonpoint contacts, namely, to the case of flat cir¬ 
cular contacts that lie on the surface of an infinite conducting 
space. This problem is of indubitable interest. Unfortunately, 
the solution that was proposed in [2] cannot be admitted cor¬ 
rect because it was found for a single flat circular contact and 
cannot be applied to a real problem, in which there are al¬ 
ways at least two current contacts of finite dimensions. As a 
result of a mutual polarization of the contacts, the boundary 
conditions are not satisfied at their borders. Therefore, the 
Schnabel solution can be used only in the case where the dis¬ 
tance between the contacts far exceeds their dimensions, and 
the approximation of point contacts then must be used. 



















anisotropic medium can be reduced to an analogous 
problem for an isotropic medium by a simple trans¬ 
formation of coordinates. 

Let there be a homogeneous conducting 
anisotropic medium that is characterized by 
the resistivity tensor p^. Let u(x\, X2, x 3 ) be the 
potential distribution in such a medium. Let us 
choose the (x 3 ,X 2 ,x 3 ) coordinate system whose 
axes are directed along the principal axes of the 
tensor p t k- In this coordinate system, the tensor 
Pik is diagonal, i.e., it has only three components 
(. Pl,P2,P3 )• 

If such a coordinate system is chosen, the equa¬ 
tion for the potential u(x\, X2, x 3 ) for an anisotropic 
medium that follows from the current-continuity 
equation (divj = 0 ) has the form 

1 d 2 u 1 d 2 u 1 d 2 u 
pi dx\ + P 2 dx 2 + P3 dx\ 

This equation differs from an ordinary Laplace 
equation, but using the linear transformation of the 
coordinate system 

Xi = f3iX* (* = 1,2,3) 

it can be easily reduced to standard Laplace equa¬ 
tion Q for an isotropic conducting medium with a 
certain resistivity p*. In this case, u(x 1 , £ 2 , # 3 ) trans¬ 
forms into u*(x\, x%, X 3 ), i.e., the potential values at 
the corresponding points of the initial anisotropic 
medium and its isotropic image are identical. Here¬ 
inafter, the asterisk in the index always designates 
the parameter of the isotropic mediumimage, which 
corresponds to the same parameter of its anisotropic 
prototype. 

The coefficients f3 1 , 02 , fi 3 , and the resistivity p* 
of the isotropic medium are chosen according to m 
so that the currents through any cross sections in 
the initial anisotropic system are equal to the cur¬ 
rents through the corresponding cross sections of its 
isotropic image, namely: 

P* = (pip2p3) 1/3 , (25) 


and 



This means that the conditional resistances mea¬ 
sured as the ratio of the potential difference between 
the potential contacts to the value of the current 
that flows between the current contacts in the initial 
anisotropic sample and its isotropic image are identi¬ 
cal, thus allowing application of the results obtained 
for isotropic samples to anisotropic samples. 


3.1 Four Collinear Contacts 


It was shown in Subsection 2.5 that for an isotropic 
sample of finite dimensions, the conditional resis¬ 
tances R\ and R 3 that are experimentally measured 
in the first and third collinear geometries are ex¬ 
pressed by the formulas 

R\ = 2 Li(a, b, d, s), R 3 = ^ L 3 (a, b, d, s), 

where the functions L\ and L 3 depend only on the 
geometrical parameters of the experiment, i.e., on 
the four dimensions a,b,d, and s, and are expressed 
by formulas (11911 and (12311 . In reality, as is seen from 
these formulas, L\ and L 3 are functions of only three 
independent arguments for which the ratios a/s, b/s, 
and a = d/s can be taken. 

For an anisotropic sample, the measured resis¬ 
tances I?i and R 3 , as was mentioned above, are the 
same as R J and R 3 for the isotropic image with the 
resistivity p* , which is determined from formula (12511 , 
and dimensions a* = (pi/p*) l / 2 a, b* = (^ 2 / P*) l ^ 2 b, 
and d* = {p 3 /p*) 1//2 d. The distance between the 
contacts in the isotropic model is s* = ( p\/p*) 1 ^ 2 s, 
if they are arranged along the x\ axis in the actual 
anisotropic sample. Hence, the expressions for R\ 
and R 3 for the anisotropic sample take the form 


Ri = ^p-L 1 (a*/s*,b*/s*,a*), 

(26) 

R3 = ^p-L 3 (a*/s*,b*/s*,a*), 


where a*/s* = a/s, b*/s* = (p 2 /pi) 1 ^ 2 b/s, and 
a* = (p 3 /piY^a, because the contacts are posi¬ 
tioned along the x\ axis. If, in addition, the condition 
Pi = P 2 (no anisotropy of p in the (. x\,X 2 ) plane is 
present) is met, the experimentally measured ratio 
of the resistances assumes the form 


R\ _ R\ 

Rs~R * 3 


L\(a/s,b/s, a*) 
L 3 (a/s,b/s,a*) 


F 13 (a/s,b/s,a*). 


(27) 

If the ratios a/s and b/s are considered as param¬ 
eters and a* as an argument, the function F 33 can 
be plotted as a function of the argument a*. As an 
example, Fig. [7Ji shows the plot of F\ 3 calculated 
from formula (12711 using formulas (11911 and (I23H for a 
sample with dimensions a = 9 mm, 6 = 3 mm, and 
s = 2.5 mm. 


Thus, by plotting the function F\ 3 {a/s,b/s,a*) = 
Fi 3 (a*) for a particular sample with the known pa¬ 
rameters a/s and b/s, we can find from this plot, 
using the measured value of the ratio R\/R 3 , the 
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3.2 Schnabel Method 




Figure 7: Plots of the functions Fi 3 (a*) = 

Li(a*)/L 3 (a*) (a) and F(a*) = G(a*)/H(a*) (b) for a 
sample with the dimensions a = 9 mm, 6 = 3 mm, and 
s = 2.5 mm. 


Measurements according to the Schnabel method 
mm should be performed in the same way as the 
procedure described above. In previous sections, ex¬ 
pressions (ED and ED were presented for the con¬ 
ditional resistances R\ and R 2 , which are measured 
for an isotropic sample in two Schnabel geometries. 
In these formulas, R\ and R 2 are expressed using 
calculated functions G ED and H ED, which de¬ 
pend only on the geometry of the experiment. As in 
Subsection 3.1, the ratio of the measured resistances 
R 1 /R 2 is a function of only three independent argu¬ 
ments: a* /s* ,6*/s*, and a* = d* /s* , where a*/s* = 
a/s, b*/s* = {p 2 /pi) 1/2 b/s, and a* = (p 3 /pi) 1/2 a, if 
the contacts are positioned along the x\ axis. 

If the condition p\ = p 2 is satisfied (no anisotropy 
of p in the (x\,X 2 ) plane is observed), the experi¬ 
mentally measured ratio of resistances assumes the 
form 


Ri 

Ro 


R/[ 

D* 

K2 


F(a/s, b/s, a*) 


G(a/s, b/s, a*) 
H(a/s, b/s, a *)' 


(29) 

Thus, by plotting the function F(a/s,b/s,a*) = 
F(a*) for a particular sample with the known pa¬ 
rameters a/s and b/s we can find the value of the 
argument a* = (p^/pi) l ^ 2 a from this plot using the 
measured ratio of the resistances R\/R 2 correspond¬ 
ing to this argument. After that, using the same 
procedure as in Subsection 3.1, the anisotropy p-z/p\ 
is found from formula (1281) . and the resistivity p\ is 
determined as 


R\d 

G(a/s,b/s,a*) H(a/s,b/s,a*) 


(30) 


value of the argument a* = (p 3 /pi) 1//2 a that cor¬ 
responds to it and, hence, determine the degree of 
anisotropy of p along the X 3 and x\ axes using the 
known geometrical parameters d and s: 


Figure [7)3 shows the plot of the function F(a*) = 
R1/R2 constructed according to formula (l29l) using 
formulae (1161) and () 12 1) for a sample with dimensions 
a = 9 mm, 6 = 3 mm, and s = 2.5 mm. 



The resistivity p\ can now be found from any for¬ 
mula ED ; 


4 Measuring the Temperature De¬ 
pendence of the Resistivity 
of Highly Oriented Pyrolitic 
Graphite 


R 1 d R%d The above-described technique was used in prac- 

Pl ~ Li(a/s, 6/s, a*) ~ L 3 (a/s,b/s,a*)’ tice to measure the resistivity of such a highly 

anisotropic material as highly oriented pyrolitic 
where the functions L\ and L 3 are calculated from graphite (HOPG) within a temperature range T = 
formulae (1191) and (1231) . and the resistivity p 3 is then 4.2-293 K using the Schnabel method. A mas- 
found from (1281) . sive HOPG sample was manufactured at the OAO 
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NIIGrafit and, according to the data of X-ray 
diffraction investigations, had an angular disorien¬ 
tation of crystallites of ~ 0.3°. In our experiment, 
the HOPG sample had the following dimensions: 
(length) a = 9 mm, (width) 6 = 3 mm, and (thick¬ 
ness) d =0.03 mm; the distance between the con¬ 
tacts was s = 2.5 mm. The plot of the function 
F(a*) = G(a*)/H(a*) in Fig. [TJo was presented, as 
an example, just for these dimensions. The sample 
was fixed in a special holder, and contacts to the sam¬ 
ple were manufactured from a 30/Um-diameter cop¬ 
per wire and attached to it with a self-solidifying 
silver paste. The experiments were performed in a 
blown through Dewar-insert for intermediate tem¬ 
peratures. The temperature stabilization and sweep 
were performed with a LakeShore model 331 temper¬ 
ature controller. 

First, the temperature dependences of the condi¬ 
tional resistances R\ and R 2 were measured. The 
obtained experimental data (Fig. [8]) for the further 
processing were represented in the form of a data ar¬ 
ray of three columns (T, R .\, R 2 ), where T is the sam¬ 
ple temperature, and i?i and R 2 are the conditional 
resistance that were measured in the first and second 
Schnabel geometries at this temperature. The num¬ 
ber of rows n in the array was determined by the 
number of experimental temperature points (in our 
case, n = 1180 at a temperature step of 0.25 K). 



Figure 8: Experimental temperature dependences of the 
conditional resistances R\ (1) and R 2 (2) in the first and 
second Schnabel geometries for a 30^m-thick HOPG sam¬ 
ple. 

The experimental data were processed in two 
stages. At the first stage, the functions G(a*), 
H(a*), and F(a*) = G(a*)/Ft (a*) were calculated 
from the known parameters (a, 6, s ) using formulas 


m and (USD. The range of values of the argument 
a* was determined by the minimum and maximum 
values of the ratio R1/R2 from the experimental data 
array (T, Ri, R 2 ). At this stage, it is necessary to 
correctly specify the depth values of the potential 
contacts £1 and £2 = 2 £i and the number of diag¬ 
onals N along which the series are summed for the 
functions G and H in order to obtain the calculation 
results with the required accuracy. At low a* val¬ 
ues, the number of diagonals N may be rather large 
(about several thousands). As the argument a* in¬ 
creases, the number of required diagonals N appre¬ 
ciably decreases. For example, the plot in Fig. [7}.) is 
constructed in the range a* = 0.5-5 with a step of 
0.25 using 19 points at £1 = O.Old and £2 = 2si. In 
this case, the number of considered diagonals along 
which summing is performed in the range a* = 0.5- 
0.8 \s N = 1600, in the range a* = 0.8-1, N = 1000, 
and at a* > 1 - N = 600. The calculation result at 
this stage is an array of numbers that contains four 
columns (a* ,G, H, F), and the number of rows in 
this array is determined by the range of values of the 
argument a* and the chosen step in this argument. 

At the second stage, the experimental data are ac¬ 
tually processed. The processing program approxi¬ 
mates the curves G(a*) and Ft (a*), which were ob¬ 
tained at the first stage, by cubic splines and con¬ 
structs the function a* = F~ l (Ri/R 2 ), which is also 
approximated by cubic splines. The processing pro¬ 
gram sequentially chooses points from the experi¬ 
mental data array (T, R\, R2), finds the ratio R1/R2, 
and then determines the value of a* = F~ 1 (R\/R 2 ) 
from the approximating formula. The anisotropy is 
determined from formula (1281) using the found value 
of a*. Subsequently, two values, p\ g and pih, are 
found for the resistivity pi from formula (1301) . into 
which the values of the functions G(a*) and Ft (a*) 
are substituted. These values are obtained using 
the above-determined approximating formulas. The 
arithmetic mean of the p\ g and p\^ values is taken 
as pi, i.e., pi = (pi g + pih)/2, and the value of p 3 is 
then found from formula (1281) . 

Figure 9 shows the results of such processing in the 
form of the temperature dependences of p\ and p 3 for 
the investigated HOPG sample and the temperature 
dependence of the value of the resistivity anisotropy 
Pz/pl- 

5 Conclusion 

Hence, both the Schnabel method and the method of 
four collinear contacts include two independent mea- 
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Figure 9: Temperature dependences of the resistivities p\ 
(1) and p 3 (2) for the HOPG sample (obtained via pro¬ 
cessing of the data presented in Fig. [HI) • The insert shows 
the temperature dependence of the resistivity anisotropy 
Pz/ Pi- 


surements of conditional resistances, which do not 
provide the determination of all three components of 
the resistivity tensor. However, if the anisotropy of a 
studied substance is such that the condition p\ = p 2 
is fulfilled, both methods allow one to find the de¬ 
gree of the anisotropy {p$/ p\) and then both values 
of pi and p 3 . The condition p\ = p -2 is satisfied 
for most substances that are interesting for practice. 
Depending on the degree of crystal anisotropy and 
its dimensions, one of these methods should be cho¬ 
sen in order to determine the desired parameters of 
a sample with the maximum accuracy. 

For example, for a sample with the dimensions 
a = 9 mm, 6 = 3 mm, and s = 2.5 mm, it fol¬ 
lows from Figs. [Tf and 7b that for “thin” samples 
(a* < 1.5), a collinear arrangement of contacts pro¬ 
vides a higher accuracy, while the Schnabel method 
should be preferred for “thick” samples (a* > 2). 

For materials in which all three resistivity compo¬ 
nents are different, other arrangements of contacts 
must be used. There are various possibilities: other 
orientations of the line of contacts, placing the po¬ 
tential contacts on the sample sides that are adjacent 
to sides on which the current contacts are positioned, 
etc. 

For samples of finite dimensions that are shaped 
as rectangular parallelepipeds, the above-described 
methods make it possible to successfully find the 
electric field potential distribution inside a sample 
and solve the formulated problem. In this case, the 


calculations are rather simple and accurate, thus al¬ 
lowing one not to use approximations in which one 
or two sample dimensions are considered infinite and 
consider an actual sample from the very beginning. 
Anyway, the technique described in this study allows 
determination of errors that arise when using the ap¬ 
proximations in which one, two, or all the three di¬ 
mensions of a sample are assumed to be infinite. 
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